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Planar quadratic differential systems occur in many areas of applied mathematics. Although 
more than one thousand papers have been written on these systems, a complete understanding 
of this family is still missing. Classical problems, and in particular, Hilbcrt's 16th problem 
[Hilbert, 1900| |Hilbert, 1902] , are still open for this family. Our goal is to make a global 
study of the family QsTiSN of all real quadratic polynomial differential systems which have a 
finite semi-elemental saddle-node and an infinite saddle-node formed by the collision of two 
infinite singular points. This family can be divided into three different subfamilies, all of them 
with the finite saddle-node in the origin of the plane with the eigenvectors on the axes and 

(A) with the infinite saddle-node in the horizontal axis, (B) with the infinite saddle-node in 
the vertical axis and (C) with the infinite saddle-node in the bisector of the first and third 
quadrants. These three subfamilies modulo the action of the affine group and time homotheties 
are three-dimensional and we give their bifurcation diagram with respect to a normal form, in 
the three-dimensional real space of the parameters of these forms. In this paper we provide 
the complete study of the geometry of the first two families, (A) and (B). The bifurcation 
diagram for the subfamily (A) yields 29 phase portraits for systems in QsnS N (A) counting 
phase portraits with and without limit cycles, while the bifurcation diagram for the subfamily 

(B) yields 16 phase portraits for systems in Q~snSN(B) under the same conditions. Case (C) 
will yield quite more cases and will have an independent paper in short. Algebraic invariants 
are used to construct the bifurcation set. The phase portraits are represented on the Poincare 
disk. The bifurcation set of Q~snSN(A) is not only algebraic due to the presence of a surface 
found numerically. All points in this surface correspond to connections of separatrices. 
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1. Introduction, brief review of the litera- 
ture and statement of results 

Here we call quadratic differential systems or simply 
quadratic systems, differential systems of the form 



x 

y 



p(x,y), 
q(x,y), 



(i) 



where p and q are polynomials over M in x and y 
such that the max(deg(p),deg((7)) = 2. To such 
a system one can always associate the quadratic 
vector field 

d d 



x = pit + 1 & ' 

ox oy 
as well as the differential equation 

qdx — pdy = 0. 



(2) 



(3) 



The class of all quadratic differential systems (or 
quadratic vector fields) will be denoted by QS. 
We can also write system ([1]) as 

x =Po +Pi(x,y) + p 2 {x,y) =p{x,y), ,^ 
V = Qo + qi(x,y) + q2(x,y) = q(x,y), 

where pi and qi are homogeneous polynomials of 
degree i in (x, y) with real coefficients withp|+g| ^ 
0. 

Even after hundreds of studies on the topology 
of real planar quadratic vector fields, it is kind of 
impossible to outline a complete characterization 
of their phase portraits, and attempting to topo- 
logically classify them, which occur rather often in 
applications, is quite a complex task. This family 
of systems depends on twelve parameters, but due 
to the action of the group G of real afhne transfor- 
mations and time homotheties, the class ultimately 
depends on five parameters, but this is still a large 
number. 

This paper is aimed at studying the class 
QsnSN of all quadratic systems possessing a fi- 
nite saddle-node srep) and an infinite saddle-node 

of type (°) SN. The finite saddle-node is a semi- 
elemental point whose neighborhood is formed 
by the union of two hyperbolic sectors and one 
parabolic sector. By a semi-elemental point we un- 
derstand a point with zero determinant of its Jaco- 
bian, but only one eigenvalue zero. These points are 
known in classical literature as semi-elementary, 
but we use the term semi-elemental introduced in 



I Artes et al., 2012] as part of a set of new definitions 
more deeply related to singular points, their multi- 
plicities and, specially, their Jacobian matrices. In 
addition, an infinite saddle-node of type (Tj SN is 
obtained by the collision of an infinite saddle with 
an infinite node. There is another type of infinite 
saddle-node denoted by (J SN which is given by 
the collision of a finite antisaddle (respectively, fi- 
nite saddle) with an infinite saddle (respectively, 
infinite node) and which will appear in some of the 
phase portraits. 

The condition of having a finite saddle-node of 
all the systems in Q~snSN implies that these sys- 
tems may have up to two other finite points. 

For a general framework of study of the 
class of all quadratic differential systems we re- 
fer to the article of Roussarie and Schlomiuk 
|Roussarie Ik Schlomiuk, 2002| . 

In this study we follow the pattern set out in 
|Artes et al, 2006] . As much as possible we shall 
try to avoid repeating technical sections which are 
the same for both papers, referring to the paper 
mentioned just above, for more complete informa- 
tion. 

In this article we give a partition of the 
classes Q'snSN(A) and Q~snSN(B). The first class 
QsnSN (A) is partitioned into 66 parts: 23 three- 
dimensional ones, 31 two-dimensional ones, 11 one- 
dimensional ones and 1 point. This partition is ob- 
tained by considering all the bifurcation surfaces 
of singularities, one related to the presence of an- 
other invariant straight line and one related to con- 
nections of separatrices, modulo "islands" (see Sec. 
[8]). The second class Q~snSN(B) is partitioned 
into 30 parts: 9 three-dimensional ones, 14 two- 
dimensional ones, 6 one-dimensional ones and 1 
point, which are all algebraic and obtained by con- 
sidering all the bifurcation surfaces. 

A graphic as defined in |Dumortier et al., 1994J 
is formed by a finite sequence of points r\, r-i, . . . ,r n 
(with possible repetitions) and non-trivial connect- 
ing orbits 7i for % = 1 , . . . , n such that 7$ has 
as a-limit set and n+i as w-limit set for i < n 
and 7 n has r n as a-limit set and r\ as w-limit set. 
Also normal orientations nj of the non-trivial or- 
bits must be coherent in the sense that if Jj—i has 
left-hand orientation then so does 7^. A polycycle 
is a graphic which has a Poincare return map. For 
more details, see |Dumortier et al, 19 94 1. 
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Theorem 1.1. There exist 29 distinct phase por- 
traits for the quadratic vector fields having a finite 
saddle-node sn(2) and an infinite saddle-node of 

type (2) SN located in the direction defined by the 
eigenvector with null eigenvalue (class QsnSN(A) ). 
All these phase portraits are shown in Figs. [7] and 
Moreover, the following statements hold: 

(a) The manifold defined by the eigenvector with 
null eigenvalue is always an invariant straight 
line under the flow; 

(b) There exist three phase portraits with limit cy- 
cles, and they are in the regions V\\, V14 and 

(c) There exist three phase portraits with graphics, 
and they are in the regions 7S\, 8S4 and I.8L1. 

Theorem 1.2. There exist 16 distinct phase por- 
traits for the quadratic vector fields having a fi- 
nite saddle-node sfi(2) and an infinite saddle-node 

of type (2) SN located in the direction defined 
by the eigenvector with non-null eigenvalue (class 
Q~snSN(B) ). All these phase portraits are shown 
in Fig. [3J Moreover, the following statements hold: 

(a) The manifold defined by the eigenvector with 
non-null eigenvalue is always an invariant 
straight line under the flow; 

(b) There exist four phase portraits with graphics, 
and they are in the regions IS4, 2S\, I.2L2 and 
1.5Li; 

(c) There exists one phase portrait with an inte- 
grable center, and it is in the region 2S±; 

(d) There exists one phase portrait with an inte- 
grable saddle, and it is in the region 2S2. 

For the class Q~snSN(A), from its 29 different 
phase portraits, 9 occur in 3-dimensional parts, 14 
in 2-dimensional parts, 5 in 1-dimensional parts 
and 1 occur in a single O-dimensional part, and for 
the class QsnSN(B), from its 16 different phase 
portraits, 5 occur in 3-dimensional parts, 7 in 2- 
dimensional parts, 3 in 1-dimensional parts and 1 
occur in a single O-dimensional part. 

In Figs, [fl [5] and E] we have denoted all the 
singular points with a small disk. We have plot- 
ted with wide curves the separatrices and we have 



added some thinner orbits to avoid confusion in 
some required cases. 

Remark 1.3. We label the phase portraits accord- 
ing to the parts of the bifurcation diagram where 
they occur. These labels could be different for two 
topologically equivalent phase portraits occurring 
in distinct parts. Some of the phase portraits in 3- 
dimensional parts also occur in some 2-dimensional 
parts bordering these 3-dimensional parts. An ex- 
ample occurs when a node turns into a focus. An 
analogous situation happens for phase portraits in 
2-dimensional (respectively, 1-dimensional) parts, 
coinciding with a phase portrait on 1-dimensional 
(respectively, O-dimensional) part situated on the 
border of it. 

The work is organized as follows. In Sec. [2] we 
describe the normal form for the families of systems 
having a finite saddle-node and an infinite saddle- 
node of type (2) SN in both horizontal and vertical 
axes. 

For the study of real planar polynomial vector 
fields two compactifications are used. In Sec. [3] we 
describe very briefly the Poincare compactification 
on the 2-dimensional sphere. 

In Sec. H]we list some very basic properties of 
general quadratic systems needed in this study. 

In Sec. [5] we mention some algebraic 
and geometric concepts that were introduced in 
ISchlomiuk et al., 200 lj |Llibre et al, 2004| involv- 
ing intersection numbers, zero-cycles, divisors, and 
T-comitants and invariants for quadratic systems 
as used by the Sibirskii school. We refer the reader 
directly to |Artes et al, 2006] where these concepts 
are widely explained. 

In Sees. [6] and [71 using algebraic invariants and 
T-comitants, we construct the bifurcation surfaces 
for the classes QJnSN(A) and Q~snSN(B), respec- 
tively. 

In Sec. [5]we comment about the possible exis- 
tence of "islands" in the bifurcation diagram. 

In Sec. [9] we introduce a global invariant de- 
noted by I, which classifies completely, up to topo- 
logical equivalence, the phase portraits we have ob- 
tained for the systems in the classes QsnSN(A) 
and QsnSN(B). Theorems 19. 81 and 19. 91 show clearly 
that they are uniquely determined (up to topologi- 
cal equivalence) by the values of the invariant I. 
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Remark 1.4- It is worth mentioning that a third 
subclass QJnSN(C) of Q'snSN must be considered. 
This subclass consists of planar quadratic systems 
with a finite saddle-node sn(2) situated as in we 
do in this work and an infinite saddle-node of type 
(2) SN in the bisector of the first and third quad- 
rants and it is currently being studied by the same 
authors. 

In I Artes et al, 1998] the authors classi- 
fied all the structurally stable quadratic pla- 
nar systems modulo limit cycles, also known as 
the codimension-zero quadratic systems (roughly 
speaking, those systems whose all singularities, fi- 
nite and infinite, are simple, with no separatrix 
connection, and where any nest of limit cycles is 
considered a single point with the stability of the 
outer limit cycle) by proving the existence of 44 
topologically different phase portraits for these sys- 
tems. The natural continuation in this idea is the 
classification of the structurally unstable quadratic 
systems of codimension-one, i.e. those systems 
which have one and only one of the following sim- 
plest structurally unstable objects: a saddle-node 
of multiplicity two (finite or infinite), a separatrix 
from one saddle point to another, and a separatrix 
forming a loop for a saddle point with its diver- 
gent non-zero. This study is already in progress 
|Artes &: Llibre, 2013| , all topological possibilities 
have already been found, some of them have already 
been proved impossible and many representatives 
have been located, but still remain some cases with- 
out candidate. One way to obtain codimension- 
one phase portraits is considering a perturbation 
of known phase portraits of quadratic systems of 
higher degree of degeneracy. This perturbation 
would decrease the codimension of the system and 
we may find a representative for a topological equiv- 
alence class in the family of the codimension-one 
systems and add it to the existing classification. 

In order to contribute to this classification, 
we study some families of quadratic systems of 
higher degree of degeneracy, e.g. systems with a 
weak focus of second order, see | Artes et al., 2006] , 
and with a finite semi-elemental triple node, see 
I Artes et al, 2013] . In this last paper, the authors 
show that, after a quadratic perturbation in the 
phase portrait Vn, the semi-elemental triple node 
is split into a node and a saddle-node and the new 



phase portrait is topologically equivalent to one of 
the topologically possible phase portrait of codi- 
mension one expected to exist. 

The present study is part of this attempt of 
classifying all the codimension-one quadratic sys- 
tems. We propose the study of a whole family of 
quadratic systems having a finite double saddle- 
node and an infinite saddle-node of type (2) SN. 
Both subfamilies reported here will not bifurcate to 
any of the codimension-one systems still missing, 
but in the subfamily QsnSN(C) will appear some 
new examples due to the fact that the complexity 
and richness of the bifurcation diagram will be the 
highest one we have already found until now. 

2. Quadratic vector fields with a finite 
saddle node Imp) and an infinite saddle 

node of type (°) SN 

A singular point r of a planar vector field X in M? 
is semi-elemental if the determinant of the matrix 
of its linear part, DX(r), is zero, but its trace is 
different from zero. 

The following result characterizes the local 
phase portrait at a semi-elemental singular point. 

Proposition 2.1. \Andronov et al., 1973\ 
Dumortier et al., 2006$ Let r = (0, 0) be an 
isolated singular point of the vector field X given 
by 

x = M(x,y), , * 

y = y + N(x,y), { 1 

where M and N are analytic in a neighborhood 
of the origin starting with at least degree 2 in the 
variables x and y. Let y = f(x) be the solu- 
tion of the equation y + N(x, y) = in a neigh- 
borhood of the point r = (0,0), and suppose that 
the function g(x) = M(x, f(x)) has the expression 
g(x) = ax a + o{x a ), where a > 2 and a / 0. So, 
when a is odd, then r = (0,0) is either an unsta- 
ble multiple node, or a multiple saddle, depending 
if a > 0, or a < 0, respectively. In the case of 
the multiple saddle, the separatrices are tangent to 
the x-axis. If a is even, the r = (0,0) is a multiple 
saddle-node, i.e. the singular point is formed by the 
union of two hyperbolic sectors with one parabolic 
sector. The stable separatrix is tangent to the posi- 
tive (respectively, negative) x-axis at r = (0,0) ac- 
cording to a < (respectively, a > 0). The two 
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unstable separatrices are tangent to the y-axis at 
r= (0,0). 

In the particular case where M and iV are real 
quadratic polynomials in the variables x and y, a 
quadratic system with a semi-elemental singular 
point at the origin can always be written into the 
form 

± = gx 2 + 2hxy + ky 2 , 

y = y + lx 2 + 2mxy + ny 2 . 

By Proposition 12.11 if g ^ 0, then we have a 
double saddle-node snp) , using the notation intro- 
duced in |Artes et ai, 2012] , 

In the normal form above, we consider the coef- 
ficient of the terms xy in both equations multiplied 
by 2 in order to make easier the calculations of the 
algebraic invariants we shall compute later. 

We note that in the normal form (|6|) we already 
have a semi-elemental point at the origin and its 
eigenvectors are (1,0) and (0, 1) which condition the 
possible positions of the infinite singular points. 

We suppose that there exists a (2) SN at some 
point at the infinity. If this point is different from 
either [1 : : 0] of the local chart Ut, or [0 : 1 : 0] 
of the local chart U2, after a reparametrization of 
the type (x,y) — > (x,ay), a € R, this point can be 
replaced at [1 : 1 : 0] of the local chart U\, that 
is, at the bisector of the first and third quadrants. 
However, if (°) SN is at [1 : : 0] or [0 : 1 : 0], 
we cannot apply this change of coordinates and it 
requires an independent study for each one of the 
cases, which in turn are not equivalent themselves 
due to the position of the infinite saddle-node with 
respect to the eigenvectors of the finite saddle-node. 

2.1. The normal form for the subclass 

Qsn~SN{A) 

The following result states the normal form for sys- 
tems in Q~snSN(A). 

Proposition 2.2. Every system with a finite 
semi-elemental double saddle-node ~sn^2) an d an in- 
finite saddle-node of typefySN located in the direc- 
tion defined by the eigenvector with null eigenvalue 
can be brought via affine transformations and time 
rescaling to the following normal form 

x = x 2 + 2hxy + ky 2 , . . 

y = y + xy + ny 2 , 



where h, k and n are real parameters. 

Proof. We start with system Q. This system al- 
ready has a finite semi-elemental double saddle- 
node at the origin (then g 7^ 0) with its eigenvec- 
tors in the direction of the axes. The first step is 
to place the point (°) SN at the origin of the local 
chart U\ with coordinates (w, z). For that, we must 
guarantee that the origin is a singularity of the flow 
in Ui, 

w = I + {—g + 2m)w + (— 2h + n)w 2 — kw 3 + wz, 
z = {—g — 2hw — kw 2 )z. 

Then, we set I = and, by analyzing the Jacobian 
of the former expression, we set m = g/2 in order to 
have the eigenvalue associated to the eigenvector on 
z = being null. Since g 7^ 0, by a time rescaling, 
we can set g = 1 and obtain the form ([?])• B 

In view that the normal form ([7]) involves the 
coefficients h, k and n, which are real, the parame- 
ter space is K 3 with coordinates (h,k,n). 

Remark 2.3. After rescaling the parameters, we 
note that system ([7]) is symmetric in relation to 
the real parameter h. Then, we shall only consider 
h > 0. 

Remark 2.4- We note that {y = 0} is an invariant 
straight line under the flow of (UJ. 

2.2. The normal form for the subclass 

QsnSN~(B) 

The following result gives the normal form for sys- 
tems in QsnSN(B). 

Proposition 2.5. Every system with a finite 
semi-elemental double saddle-node sn(2) and an in- 
finite saddle-node of typefySN located in the direc- 
tion defined by the eigenvector with non-null eigen- 
value can be brought via affine transformations and 
time rescaling to the following normal form 

x = x 2 + 2hxy, . . 

y = y + lx 2 + 2mxy + 2hy 2 , 

where h, I and m are real parameters. 

Proof. Analogously to Proposition 12.21 we start 
with system (JH|) , but now we want to place the point 
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(ASN at the origin of the local chart XJ%. By fol- 
lowing the same steps, we set k = 0, n = h/2, g = 1 
and we obtain the form (JSj) - B 

In view that the normal form ([7|) involves the 
coefficients h, I and m, which are real, the parame- 
ter space is M 3 with coordinates (h,l,m). 

Remark 2.6. After rescaling the parameters, we 
note that system (|8|) is symmetric in relation to 
the real parameter h. Then, we will only consider 
h > 0. 

Remark 2.7. We note that {x = 0} is an invariant 
straight line under the flow of ([8]). 

3. The Poincare compactification and the 
complex (real) foliation with singulari- 
ties on CP 2 (MP 2 ) 

A real planar polynomial vector field £ can be com- 
pactified on the sphere as follows. Consider the 
x, y plane as being the plane Z = 1 in the space 
M 3 with coordinates X, Y, Z. The central pro- 
jection of the vector field £ on the sphere of ra- 
dius one yields a diffeomorphic vector field on the 
upper hemisphere and also another vector field on 
the lower hemisphere. There exists (for a proof see 
|Gonzales, 1969| ) an analytic vector field cp(£) on 
the whole sphere such that its restriction on the 
upper hemisphere has the same phase curves as the 
one constructed above from the polynomial vector 
field. The projection of the closed northern hemi- 
sphere H + of S 2 on Z = under (X, Y, Z) -> 
(X, Y) is called the Poincare disc. A singular point 
q of cp(£) is called an infinite (respectively, finite) 
singular point if q € S , the equator (respectively, 
q € S 2 \S 1 ). By the Poincare compactification of 
a polynomial vector field we mean the vector field 
cp(£) restricted to the upper hemisphere completed 
with the equator. 

Ideas in the remaining part of this section 
go back to Darboux's work |Darboux, 1878] . Let 
p(x,y) and q(x,y) be polynomials with real coeffi- 
cients. For the vector field 



or equivalently for the differential system 

x = p{x,y), y = q{x,y), (10) 



we consider the associated differential 1-form 
01 = q(x, y)dx—p(x, y)dy, and the differential equa- 
tion 

01=0. (11) 

Clearly, equation (jlip defines a foliation with sin- 
gularities on C 2 . The affine plane C 2 is com- 
pactified on the complex projective space CP 2 = 
(C 3 \ {0})/ ~, where (X, Y, Z) ~ (X', Y', Z') if and 
only if {X, Y,Z) = \{X' ,Y' , Z') for some complex 
A / 0. The equivalence class of (X, Y, Z) will be 
denoted by [X : Y : Z]. 

The foliation with singularities defined by equa- 
tion (|lip on C 2 can be extended to a foliation with 
singularities on CP 2 and the 1-form 0i can be ex- 
tended to a meromorphic 1-form on CP 2 which 
yields an equation = 0, i.e. 

A(X,Y, Z)dX + B(X,Y, Z)dY + C{X,Y, Z)dZ = 0, 

(12) 

whose coefficients A, B, C are homogeneous poly- 
nomials of the same degree and satisfy the relation: 

A(X, Y, Z)X + B(X, Y, Z)Y + C(X, Y, Z)Z = 0, 

(13) 

Indeed, consider the map i : C 3 \ {Z = 0} — > C 2 , 
given by i(X, Y, Z) = (X/Z, Y/Z) = (x, y) and sup- 
pose that max{deg(p), deg(g)} = m > 0. Since 
x = X/Z and y = Y/Z we have: 

dx = (ZdX - XdZ) /Z 2 , dy = (ZdY - YdZ) /Z 2 , 

the pull-back form i*{0\) has poles at Z = and 
yields the equation 

i*(0i) =q(X/Z,Y/Z)(ZdX - XdZ)/Z 2 

-p(X/Z,Y/Z)(ZdY - YdZ)/Z 2 = 0. 

Then, the 1-form = Z m+2 i*(0 1 ) in C 3 \ {Z ^ 0} 
has homogeneous polynomial coefficients of degree 
m + 1, and for Z = the equations = and 
i*(0i) = have the same solutions. Therefore, the 
differential equation = can be written as (|12|) . 
where 

A(X, Y, Z) =ZQ(X, Y, Z) = Z m+1 q(X/Z, Y/Z), 
B(X, Y,Z) = - ZP(X, Y, Z) = -Z m+1 p(X/Z, Y/Z), 
C(X, Y, Z) =YP(X, Y, Z) - XQ(X, Y, Z). 

(14) 

Clearly A, B and C are homogeneous polyno- 
mials of degree m + 1 satisfying (|13D . 
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In particular, for our quadratic systems ([7]), A, 
B and C take the following forms 

A(X, Y, Z) =YZ(X -nY + Z) 

B(X, Y,Z) = - (X 2 + 2hXY + kY 2 )Z, 

C(X, Y, Z) =Y(2hXY - nXY + kY 2 - XZ). 

(15) 

and for our quadratic systems A, B and C take 
the following forms 

A(X, Y, Z) =Z{IX 2 + mXY + 2hY 2 + YZ), 

B(X, Y,Z)=- X(X + 2hY)Z, 

C(X, Y, Z) =X - X(IX 2 -XY + 2mXY + YZ). 

(16) 

We note that the straight line Z = is always 
an algebraic invariant curve of this foliation and 
that its singular points are the solutions of the sys- 
tem: A(X,Y,Z) = B(X,Y,Z) = C(X,Y,Z) = 0. 
We note also that C(X,Y, Z) does not depend on 
b. 

To study the foliation with singularities defined 
by the differential equation (fl~2j) subject to (fl"3j) 
with A, B, C satisfying the above conditions in the 
neighborhood of the line Z = 0, we consider the 
two charts of CP 2 : (u,z) = (Y/X,Z/X), I/O, 
and (v,w) = (X/Y, Z/Y), Y ^ 0, covering this 
line. We note that in the intersection of the charts 
(x, y) = (X/Z,Y/Z) and (u, z) (respectively, (v,w)) 
we have the change of coordinates x = 1/z, y = u/z 
(respectively, x = v/w, y = 1/w). Except for the 
point [0 : 1 : 0] or the point [1 : : 0], the foliation 
defined by equations (fl~2|) . (fl"3|) with A, B, C as in 
(|14p yields in the neighborhood of the line Z = 
the foliations associated with the systems 

u =uP(l,u, z) - Q(l, u, z) = (7(1, u, z), ^ 
z =zP(l, u, z), 

or 

v =vQ(v,l,w) - P(v,l,w) = -C(v,l,w), 
w =wP(v, 1, W). 

In a similar way we can associate a real foliation 
with singularities on MP 2 to a real planar polyno- 
mial vector field. 

4. A few basic properties of quadratic sys- 
tems relevant for this study 

We list below results which play a role in the 
study of the global phase portraits of the real pla- 



nar quadratic systems ([T]) having a semi-elemental 
triple node. 

The following results hold for any quadratic 
system: 

(i) A straight line either has at most two (finite) 
contact points with a quadratic system (which 
include the singular points), or it is formed by 
trajectories of the system; see Lemma 11.1 of 
|Ye et al, 1986| . We recall that by definition 
a contact point of a straight line L is a point of 
L where the vector field has the same direction 
as L, or it is zero. 

(ii) If a straight line passing through two real fi- 
nite singular points r\ and T2 of a quadratic 
system is not formed by trajectories, then it is 
divided by these two singular points in three 
segments oofT, TTT2 and r^oo such that the 
trajectories cross oorT and r200 in one direc- 
tion, and they cross in the opposite di- 
rection; see Lemma 11.4 of |Ye et al., 1986| . 

(hi) If a quadratic system has a limit cycle, then 
it surrounds a unique singular point, and this 
point is a focus; see |Coppel, 1966 . 

(iv) A quadratic system with an invariant straight 
line has at most one limit cycle; see 
|Coll Llibre, 1988] . 

(v) A quadratic system with more than one in- 
variant straight line has no limit cycle; see 
IBautin, 1954| . 

Proposition 4.1. Any graphic or degenerate 
graphic in a real planar polynomial differential 
system must either 

1 ) surround a singular point of index greater than 
or equal to +1, or 

2) contain a singular point having an elliptic sector 
situated in the region delimited by the graphic, or 

3) contain an infinite number of singular points. 
Proof. See the proof in |Artes et al., 1998] . ■ 

5. Some algebraic and geometric concepts 

In this article we use the concept of intersection 
number for curves (see |Fulton, 1969| ). For a quick 
summary see Sec. 5 of |Artes et al., 2006] . 
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We shall also use the concepts of zero- 
cycle and divisor (see |Hartshorne, 1977] ) 
as specified for quadratic vector fields in 
|Schlomiuk et al, 2001| . For a quick summary see 



Sec. 6 of lArtes et al, 2006 



We shall also use the concepts of algebraic in- 
variant and T-comitant as used by the Sibirskii 
school for differential equations. For a quick sum- 
mary see Sec. 7 of | Artes et al., 2006] . 

In the next two sections we describe the alge- 
braic invariants and T-comitants which are relevant 
in the study of families (0), see Sec. El and ( 
Sec. El 



sec 



6. The bifurcation diagram of the systems 
in Qsn~SN(A) 

We recall that, in view that the normal form dTJ) 
involves the coefficients h, k and n, which are real, 
the parameter space is M 3 with coordinates (h, k, n). 

6.1. Bifurcation surfaces due to the 
changes in the nature of singularities 

For systems ([7]) we will always have (0, 0) as a finite 
singular point, a double saddle-node. 

From Sec. 7 of |Artes et al, 2008| we get the 
formulas which give the bifurcation surfaces of sin- 
gularities in M 12 , produced by changes that may oc- 
cur in the local nature of finite singularities. From 



Schlomiuk et al., 2005 we get equivalent formulas 
for the infinite singular points. These bifurcation 
surfaces are all algebraic and they are the follow- 
ing: 



Bifurcation surfaces in 
ties of singularities 



due to multiplici- 



(5i) This is the bifurcation surface due to mul- 
tiplicity of infinite singularities as detected by 
the coefficients of the divisor D^{P,Q;Z) = 
T,wa{z =0}nCP 2j w{P,Q)W, (here I W {P,Q) de- 
notes the intersection multiplicity of P = with 
Q = at the point W situated on the line at infin- 
ity, i.e. Z = 0) whenever deg({D R (P,Q; Z))) > 0. 
This occurs when at least one finite singular point 
collides with at least one infinite point. More 
precisely this happens whenever the homogenous 
polynomials of degree two, p2 and q2, in p and 
q have a common root. In other words whenever 



[i = ReSa;(p2) ?2)/y 4 = 0. The equation of this sur- 
face is 

fi = k — 2hn + n 2 = 0. 

Since this family already has a saddle- 
node at the origin, the invariant B, as defined in 
|Artes et al, 200"6] , is always zero. Since our nor- 
mal form does not allow the existence of a finite 
singular point of multiplicity 3, the only possible 
bifurcation related to collisions of finite singulari- 
ties with themselves is whether the other two finite 
singularities are either real, or complex, or form a 
double one. This phenomenon is captured by the 
T-comitant T as proved in |Artes et al, 200*8] . The 
equation of this surface is 

T = h 2 - k = 0. 

(55) Since this family already has a saddle- 
node at infinity, the invariant r/, as defined in 
|Artes et al, 2006] , is always zero. In this sense, 
we have to consider a bifurcation related to the 
existence of either the double infinite singularity 
(?i)SN plus a simple one, or a triple one. This phe- 
nomenon is ruled by the T-comitant M as proved 
in |Schlomiuk et al., 20051 [Artes et al, "2012] . The 
equation of this surface is 

M = 2h - n = 0. 

The surface of C°° bifurcation points due to 
a strong saddle or a strong focus changing 
the sign of their traces (weak saddle or weak 
focus) 

(1S2) This is the bifurcation surface due to the weak- 
ness of finite singularities, which occurs when the 
trace of a finite singular point is zero. The equation 
of this surface is given by 



-Ah 2 + 4k + n 2 = 0, 



where Ta is defined in |Vulpe, 2011| . This Ta is an 
invariant. 

This bifurcation can produce a topological 
change if the weak point is a focus or just a C°° 



The numbers attached to these bifurcations surfaces do 
not appear here in increasing order. We just kept the same 
enumeration used in | Artes et al., 2006] to maintain coher- 
ence even though some of the numbers in that enumeration 
do not occur here. 
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change if it is a saddle, except when this bifurcation 
coincides with a loop bifurcation associated with 
the same saddle, in which case, the change may 
also be topological. 

The surface of C°° bifurcation due to a node 
becoming a focus 

(Sq) This surface will contain the points of the pa- 
rameter space where a finite node of the system 
turns into a focus. This surface is a C°° but not 
a topological bifurcation surface. In fact, when 
we only cross the surface (Sq) in the bifurcation 
diagram, the topological phase portraits do not 
change. However, this surface is relevant for iso- 
lating the regions where a limit cycle surrounding 
an antisaddle cannot exist. Using the results of 
|Artes et al, 2008] , the equation of this surface is 
given by W4 = 0, where 

W 4 = - 48/i 4 + 32h 2 k + 16k 2 + 64/i 3 n 
- QAhkn - 2Ah 2 n 2 + 2Akn 2 + n A . 

Bifurcation surface in ]R 3 due to the presence 
of another invariant straight line 

(Ss) This surface will contain the points of the pa- 
rameter space where another invariant straight line 
appears apart from {y = 0}. This surface is split 
in some regions. Depending on these regions, the 
straight line may contain connections of separatri- 
ces from different saddles or not. So, in some cases, 
it may imply a topological bifurcation and, in oth- 
ers, just a C°° bifurcation. The equation of this 
surface is given by 

Het = h = 0. 

These, except (Ss), are all the bifurcation sur- 
faces of singularities of systems (|7J) in the parame- 
ter space and they are all algebraic. We shall dis- 
cover another bifurcation surface not necessarily al- 
gebraic and on which the systems have global con- 
nection of separatrices different from that given by 
(5s). The equation of this bifurcation surface can 
only be determined approximately by means of nu- 
merical tools. Using arguments of continuity in the 
phase portraits we can prove the existence of this 
not necessarily algebraic component in the region 
where it appears, and we can check it numerically. 
We will name it the surface (57). 



Remark 6.1. Even though we can draw a 3- 
dimensional picture of the algebraic bifurcation sur- 
faces of singularities in M 3 , it is pointless to try to 
see a single 3-dimensional image of all these bifur- 
cation surfaces together in the space M 3 . As we 
shall see later, the full partition of the parameter 
space obtained from all these bifurcation surfaces 
has 66 parts. 

By the previous remark we shall foliate the 3- 
dimensional bifurcation diagram in R 3 by planes 
h = ho, ho constant. We shall give pictures of the 
resulting bifurcation diagram on these planar sec- 
tions on an affine chart on M 2 . In order to detect the 
key values for this foliation, we must find the values 
of parameters where the surfaces have singularities 
and/or intersect to each other. As we mentioned 
before, we will be only interested in non-negative 
values of h to construct the bifurcation diagram. 

The following set of sixteen results study the 
singularities of each surface and the simultaneous 
intersection points of the bifurcation surfaces, or 
the points or curves where two bifurcation surfaces 
are tangent. 

As the final bifurcation diagram is quite com- 
plex, it is useful to introduce colors which will be 
used to talk about the bifurcation points: 

(a) the curve obtained from the surface (Si) is 
drawn in blue (a finite singular point collides 
with an infinite one); 

(b) the curve obtained from the surface (S2) is 
drawn in yellow (when the trace of a singular 
point becomes zero); 

(c) the curve obtained from the surface (S3) is 
drawn in green (two finite singular points col- 
lide); 

(d) the curve obtained from the surface (S5) is 
drawn in red (three infinite singular points col- 
lide); 

(e) the curve obtained from the surface (Sq) is 
drawn in black (an antisaddle is on the edge of 
turning from a node to a focus or vice versa); 

(f) the curve obtained from the surface (1S7) is 
drawn in purple (the connection of separatri- 
ces); and 
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(g) the curve obtained from the surface (Ss) is 
also drawn in purple (presence of an invariant 
straight line) . We draw it as a continuous curve 
if it implies a topological change or as a dashed 
curve if not. 

We use the same color for (57) and (Ss) since 
both surfaces deal with connections of separatrices 
mostly. 

Lemma 6.2. Concerning the singularities of the 
surfaces, it follows that: 

(i) (Si), (S2), (S3), (S5) and (Ss) have no sin- 
gularities; 

(ii) (Sq) has singularities at the point (0,0,0) and 
along the straight line (h,0,2h). 

Proof. It is easy to see that the gradient vectors 
of each one of the surfaces (S\), (S2), (S3), (S5) 
and (Ss) are never null for all (h,k,n) G R 3 ; so 
(i) is proved. In order to prove (ii) we compute the 
gradient of W4 and we verify that it is null whenever 
h = k = n = and along the straight line (h, 0, 2h). 

U 

Lemma 6.3. The surfaces (S\) and (S2) intersect 
along the straight line (h, 0, 2h) and along the curve 
(h,8h 2 /9,2h/3), for all heR. 

Proof. By solving simultaneously both equations 
of the surfaces (<Si) and (52) for all h 7^ 0, we 
obtain the straight line (h, 0, 2h) and the curve 
(h, 8h 2 /9, 2/i/3); if h = 0, the only simultaneous 
solution is the origin. ■ 

The proofs of the next lemmas are analogous 
to Lemma 16.31 except when a different proof is in- 
cluded. 

Lemma 6.4. If h = 0, the surfaces (S\) and (S3) 
intersect at the origin and, if h 7^ 0, both surfaces 
intersect along the curve (n,n 2 ,n), n G R. 

Lemma 6.5. For all h G R, the surfaces (S\) and 
(S§) intersect along the straight line (h,0,2h). 

Lemma 6.6. For all h G R, the surfaces (S\) and 
(Sq) intersect along the straight line (h,0,2h) and 
along the curve (/i, 48/i 2 /49, 6h/7). 



Lemma 6.7. For all h G R, the surfaces (S\) and 
(Ss) intersect along the straight line (h,0,2h). 

Lemma 6.8. For all h G R, the surfaces (S2) and 
(S3) intersect along the curve (h,h 2 ,0). 

Lemma 6.9. For all h G R, the surfaces (S2) and 
(S5) intersect along the straight line (h,0,2h). 

Lemma 6.10. For all the surfaces (S2) and 

(Sq) intersect along the straight line (h,0,2h) and 
along the curve (h,h 2 ,0). 

Lemma 6.11. For all h G R, the surfaces (S2) and 
(S$) intersect along the straight lines (h,0,—2h) 
and (h, 0, 2h). 

Lemma 6.12. For all h G R, the surfaces (S3) and 
(S5) intersect along the curve (h,h 2 ,2h). 

Lemma 6.13. For all h G R, the surfaces (S3) and 
(Sq) intersect along the curve (h,h 2 ,0). 

Lemma 6.14. Ifh = 0, the surfaces (S3) and (S&) 
intersect along the straight line (0, 0, n), n G R, 
and, if h 7^ 0, they have no intersection. 

Proof. By restricting the equations of both surfaces 
to h = and solving them simultaneously we obtain 
the straight line (0, 0,n), n G R. For all h 7^ 0, the 
equation have no simultaneous solutions. ■ 

Lemma 6.15. For all the surfaces (S5) and 

(Sq) intersect along the straight line (h,0,2h). 

Lemma 6.16. For all ft€R, the surfaces (S5) and 
(Ss) intersect along the straight line (h,0,2h). 

Lemma 6.17. For all h G R, the surfaces (Sq) and 
(Ss) intersect along the straight lines (h,0,—6h) 
and (h,0,2h). 

Now we shall study the bifurcation diagram 
having as reference the values of h where significant 
phenomena occur in the behavior of the bifurcation 
surfaces. As there is not any other critical value of 
h, except h = 0, this is the only value where the be- 
havior of the bifurcation surfaces changes critically. 
Recalling we are considering only non-negative val- 
ues of h, we shall choose a positive value to be a 
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Fig. 4. Slice of the parameter space for ([7]) when h = 0. 

generic case. 

We take then the values: 



The value ho corresponds to an explicit value 
of h for which there is a bifurcation in the behavior 
of the systems on the slices. The value hi is just 
an intermediate point we call by a generic value of 
h (see Figs. [3] and [5]). 

We now describe the labels used for each part. 
The subsets of dimensions 3, 2, 1 and 0, of the parti- 
tion of the parameter space will be denoted respec- 
tively by V, 5, L and P for Volume, Surface, Line 
and Point, respectively. The surfaces are named 
using a number which corresponds to each bifurca- 
tion surface which is placed on the left side of the 
letter 5. To describe the portion of the surface we 
place an index. The curves that are intersection 
of surfaces are named by using their corresponding 
numbers on the left side of the letter L, separated 
by a point. To describe the segment of the curve 
we place an index. Volumes and Points are sim- 
ply indexed (since three or more surfaces may be 
involved in such an intersection). 

We consider an example: the surface (Si) splits 
into 6 different two-dimensional parts labeled from 
lSi to 156, plus some one-dimensional arcs labeled 
as l.iLj (where i denotes the other surface inter- 
sected by (S\) and j is a number), and some zero- 
dimensional parts. In order to simplify the labels 
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Fig. 5. Slice of the parameter space for Q when h = l. 

in Figs. [8] and [9] we see VI which stands for the 
TpjX notation Vi. Analogously, 1S1 (respectively, 
1.2L1) stands for 15i (respectively, 1.2Li). And 
the same happens with other pictures. 

All the bifurcation surfaces intersect on h = 
0. In fact, the equations of surfaces (Si) and (52 ) 
are the parabolas k + n 2 = and 4k + n 2 = 0, 
respectively, when restricted to the plane h = 0; 
the equations of (53), (S$) and (5s), restricted to 
h = 0, are the straight lines k = 0, n = and k = 0, 
respectively; and surface (56 ) is the quartic 16k 2 + 
24kn 2 + n 4 = whose picture is the union of two 
parabolas with intersection at the origin. Finally, 
we note that all the elements above intersect at the 
origin of the bifurcation diagram when h = 0. 

As an exact drawing of the curves produced by 
intersecting the surfaces with slices gives us very 
small regions which are difficult to distinguish, 
and points of tangency are almost impossible 
to recognize, we have produced topologically 
equivalent pictures where regions are enlarged 
and tangencies are easy to observe. The reader 
may find the exact pictures in the web page 
http:/ /mat .uab.es /^artes /articles / qvfsn2SN02 / 
qvfsn2SN02.html. 

If we consider the value h = 1 , some changes in 
the bifurcation diagram happen. On the one hand, 
all the surfaces preserve their geometrical behav- 
ior, that is, surfaces (Si) and (52) remain parabolas 
(k — 2n + n 2 = and —4 + 4kn 2 = 0, respectively); 
surfaces (53), (S5) and (5s) remain straight lines 
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(k — 1 = 0, n — 2 = and k = 0, respectively); 
and surface (Sq) remains a quartic (with equation 
-48 + 32/c + 1 6 A; 2 + 64n - 64/cn - 24n 2 + 24 kn 2 + n 4 = 
0) whose picture is the union of two curves with in- 
tersection at the point (1,0, 2). On the other hand, 
compared to the case when h = 0, there exist more 
intersection points among the surfaces and a new 
region appears between (S3) and (5s). All other re- 
gions, except this new "middle" one, remain topo- 
logically equivalent to the regions present in the 
case when h = 0. 

We recall that the black surface (Sq) (or W4) 
means the turning of a finite antisaddle from a node 
to a focus. Then, according to the general results 
about quadratic systems, we could have limit cycles 
around such point. 

Remark 6.18. Wherever two parts of equal dimen- 
sion d are separated only by a part of dimension 
d — 1 of the black bifurcation surface (Sq), their 
respective phase portraits are topologically equiva- 
lent since the only difference between them is that 
a finite antisaddle has turned into a focus without 
change of stability and without appearance of limit 
cycles. We denote such parts with different labels, 
but we do not give specific phase portraits in pic- 
tures attached to Theorems 1 1 . 1 1 and 1 1 . 2 1 for the parts 
with the focus. We only give portraits for the parts 
with nodes, except in the case of existence of a limit 
cycle or a graphic where the singular point inside 
them is portrayed as a focus. Neither do we give 
specific invariant description in Sec. [9] distinguish- 
ing between these nodes and foci. 

6.2. Bifurcation surfaces due to connec- 
tions 

We now describe for each set of the partition on 
h = 1 the local behavior of the flow around all the 
singular points. Given a concrete value of param- 
eters of each one of the sets in this slice we com- 
pute the global phase portrait with the numerical 
program P4 |Dumortier et al, 2006| . It is worth 
mentioning that many (but not all) of the phase 
portraits in this paper can be obtained not only 
numerically but also by means of perturbations of 
the systems of codimension one higher. 

In this slice we have a partition in 2- 
dimensional regions bordered by curved polygons, 
some of them bounded, others bordered by infinity. 



Provisionally, we use low-case letters to describe 
the sets found algebraically so as not to interfere 
with the final partition described with capital let- 
ters. For each 2-dimensional region we obtain a 
phase portrait which is coherent with those of all 
their borders, except in one region. Consider the 
set v\ in Fig. [5j In it we have only a saddle-node 
as finite singularity. When reaching the set 2.3l\, we 
are on surfaces (1S2), (S3) and (Sq) at the same time; 
this implies the presence of one more finite singular- 
ity (in fact, it is a cusp point) which is on the edge 
of splitting itself and give birth to finite saddle and 
antisaddle. Now, we consider the segments 2s2 and 
2s 3 . By the Main Theorem of |Vulpe, 2011] , the 
corresponding phase portraits of these sets have a 
first-order weak saddle and a first-order weak fo- 
cus, respectively. So, on 2S3 we have a Hopf bifur- 
cation. This means that either in V5 or v\q we must 
have a limit cycle. In fact, it is in V5. On the other 
hand, as we have a weak saddle on 2s2 and it is 
not detected a bifurcation surface intersecting this 
subset of loop type, neither its presence is forced 
to keep the coherence, its corresponding phase por- 
trait is topologically equivalent to the portraits of 
t>4 and v$. Since in v§ we have a phase portrait topo- 
logically equivalent to the one on 2s2 (without limit 
cycles) and a phase portrait with limit cycles, this 
region must be split into two regions separated by 
a new surface (S7) having at least one element 7 Si 
such that one region has limit cycle and the other 
does not, and the border 7S± must correspond to 
a connection between separatrices. After numerical 
computations we check that it is the region v§ the 
one which splits into V5 without limit cycles and 
V\\ with one limit cycle. 

The next result assures us the existence of limit 
cycle in any representative of the subset v\± and it 
is needed to complete the study of 7Si . 

Lemma 6.19. In v\4 there is always one limit cy- 
cle. 

Proof. We see that the subset U14 is characterized 
by a < 0, 74 < 0, W4 < 0, M > 0, T > 0, k > 
and n > 0. Any representative of v 14 has the finite 
saddle-node at the origin with its eigenvectors on 
the axes and two more finite singularities, a focus 
and a node (the focus is due to W4 < 0). We claim 
that these two other singularities are placed in sym- 
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Fig. 6. The local behavior around each of the finite and 
infinite singularities of any representant of v\± . The red 
arrow shows the sense of the flow along the y-axis and 
the blue points are the focus and the node with same 
stability. 

metrical quadrants with relation to the origin (see 
Fig. In fact, by computing the exact expression 
of each singular point (x\, y±) and (x2, J/2) and mul- 
tiplying their x -coordinates and y-coordinates we 
obtain k/fj, and 1/fx, respectively, which are always 
negative since k > and fi < in v\±. Besides, each 
one of them is placed in an even quadrant since the 
product of the coordinates of each antisaddles is 
never null and any representative gives a negative 
product. Moreover, both antisaddles have the same 
stability since the product of their traces is given by 
li/Ta which is always positive in vi^. 

The infinite singularities of systems in v\a are 
the saddle-node (^\SN (recall the normal form (|5])) 
and a saddle. In fact, the expression of the sin- 
gular points in the local chart U\ are (0, 0) and 
((— 2h + n)/k,0). We note that the determinant of 
the Jacobian matrix of the flow in U\ at the second 
singularity is given by 

(2h - n) 2 (2hn - k - n 2 ) M 2 li 

which is negative since in < in V14. Besides, this 
pair of saddles are in the second and the fourth 
quadrant because its first coordinate (—2h + n)/k = 
—M/k is negative since M > and k > in 1*14. 

We also note that the flow along the y-axis is 
such that x > 0. 

Since we have a pair of saddle points in the even 
quadrants, each of the finite antisaddles is in an 
even quadrant, no orbit can enter into the second 
quadrant and no orbit may leave the fourth one 



and, in addition, these antisaddles, a focus and a 
node, have the same stability, any phase portrait 
in v 14 must have at least one limit cycle in any of 
the even quadrants. Moreover, the limit cycle is 
in the second quadrant, because the focus is there 
since a saddle-node is born in that quadrant at 3s\, 
splits in two points when entering v% (both remain 
in the same quadrant since X1X2 = k/fj, < and 
yiU2 = < 0), the node turns into focus at 6S2 
and the saddle moves to infinity at 1^2 appearing 
as node at the fourth quadrant when entering V14. 
Furthermore, by the statement (jrv|) of Sec. 01 it 
follows the uniqueness of the limit cycle in vu- ■ 

Now, the following result states that the seg- 
ment which splits the subset v$ into the regions V5 
and V11 has its endpoints well-determined. 

Proposition 6.20. The endpoints of the surface 
7S\ are 2.3l\, intersection of surfaces (S2) and 
(S3), and I.8/1, intersection of surfaces (S\) and 
(Ss). 

Even though the next proof will be done in the 
concrete case when h = 1, it can be easily extended 
for the generical case of surfaces and curves. We 
can visualize the image of this surface in the plane 
h = lm Fig. El 

Proof We consider h = 1 and we write r\ = (1, 0, 2) 
and r2 = (1, 1,0) for 2.3l\ and I.8Z1 , respectively. If 
the starting point were any point of the segments 
2s2 and 2S3, we would have the following incoher- 
ences: firstly, if the starting point of 7S\ were on 
2s2, a portion of this subset must refer to a Hopf 
bifurcation since we have a limit cycle in Vu; and 
secondly, if this starting point were on 2ss, a por- 
tion of this subset must not refer to a Hopf bifurca- 
tion which contradicts the fact that on 2^3 we have 
a first-order weak focus. Finally, the ending point 
must be r2 because, if it were located on 8S3, we 
would have a segment between this point and I.8/1 
along surface (Ss) with two invariant straight lines 
and one limit cycle, which contradicts the statement 
(jvj) of Sec. HI and if it were on 1^2, we would have a 
segment between this point and I.8/1 along surface 
(<Si) without limit cycle which is not compatible 
with Lemma 16.191 since ll = does not produce a 
graphic. ■ 
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We show the sequence of phase portraits along 
these subsets in Fig. [7J 

We cannot be totally sure that this is the 
unique additional bifurcation curve in this slice. 
There could exist others which are closed curves 
which are small enough to escape our numerical re- 
search, but the located one is enough to maintain 
the coherence of the bifurcation diagram. We re- 
call that this kind of studies are always done mod- 
ulo "islands". For all other two-dimensional parts 
of the partition of this slice whenever we join two 
points which are close to two different borders of 
the part, the two phase portraits are topologically 
equivalent. So we do not encounter more situations 
than the one mentioned above. 

As we vary h in (0, oo), the numerical research 
shows us the existence of the phenomenon just de- 
scribed, but for h = 0, we have not found the same 
behavior. 

In Figs. [8] and [9] we show the complete bifurca- 
tion diagrams. In Sec. [U]the reader can look for the 
topological equivalences among the phase portraits 
appearing in the various parts and the selected no- 
tation for their representatives in Figs. [1] and 

7. The bifurcation diagram of the systems 
in QsnSN(B) 

We recall that, in view that the normal form (JSj) 
involves the coefficients h, I and m, which are real, 
the parameter space is M 3 with coordinates (h, I, m). 

Before we describe all the bifurcation surfaces 
for Q~snSN(B), we prove the following result which 
gives conditions on the parameters for the presence 
of either a finite star node n* (whenever any two dis- 
tinct non-trivial integral curves arrive at the node 
with distinct slopes), or a finite dicritical node n d (a 
node with identical eigenvalues but Jacobian non- 
diagonal). 

Lemma 7.1. Systems always have an*, ifm = 
and h ^ 0, or a n d , otherwise. 

Proof. We note that the singular point (0, —l/2h) 
has its Jacobian matrix given by 

( -1 \ 
V -m/h -1 J ' 



7.1. Bifurcation surfaces due to the 
changes in the nature of singularities 

For systems ([8]) we will always have (0, 0) as a fi- 
nite singular point, a double saddle-node. Besides, 
the needed invariants here are the same as in the 
previous system except surfaces (Sq) and (Ss); so 
that we shall only give the geometrical meaning and 
their equations plus a deeper discussion on surface 
(Sq). For further information about them, see Sec. 

El 
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Fig. 7. Sequence of phase portraits in slice h = 1 from v\ to I.8Z1. We start from v\. When crossing 2.3/i, we may 
choose at least seven "destinations": 6S2, Vi, 2s2, ^5, 2ss, Vio and 6S3. In each one of these subsets, but 115, we 
obtain only one phase portrait. In v§ we find (at least) three different ones, which means that this subset must be 
split into (at least) three different regions whose phase portraits are V5, 7 Si and V\i- And then we shall follow the 
arrows to reach the subset I.8/1 whose corresponding phase portrait is I.8L1. 
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Fig. 9. Complete bifurcation diagram of QsnSN(A) for slice h=l. 
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Bifurcation surfaces in IR 3 due to multiplici- 
ties of singularities 

(Si) This is the bifurcation surface due to multi- 
plicity of infinite singularities. This occurs when at 
least one finite singular point collides with at least 
one infinite point. The equation of this surface is 

H = 4/i 2 (l + 2hl - 2m) = 0. 

(^3) This is the bifurcation surface due to the col- 
lision of the other two finite singularities different 
from the saddle-node. The equation of this surface 
is given by 

T = —h 2 = 0. 

It only has substantial importance when we con- 
sider the plane h = 0. 

(55) This is the bifurcation surface due to the col- 
lision of infinite singularities, i.e. when all three in- 
finite singular points collide. The equation of this 
surface is 

M= (-1 + 2m) 2 = 0. 

The surface of C°° bifurcation due to a strong 
saddle or a strong focus changing the sign of 
their traces (weak saddle or weak focus) 

(£2) This is the bifurcation surface due to the weak- 
ness of finite singularities, which occurs when their 
trace is zero. The equation of this surface is given 
by 

% = -I6/1 3 / = 0. 

The surface of C°° bifurcation due to a node 
becoming a focus 

(Sq) Since W4 is identically zero for all the bifurca- 
tion space, the invariant that captures if a second 
point may be on the edge of changing from node to 
focus is W3. The equation of this surface is given 
by 

W 3 = 64/i 4 (l + 2hl + h 2 l 2 - 2m) = 0. 

These are all the bifurcation surfaces of singu- 
larities of the systems (JHJ) in the parameter space 
and they are all algebraic. We are not meant to 
discover any other bifurcation surface (neither non- 
algebraic nor algebraic one) due to the fact that in 
all the transitions we make among the parts of the 



bifurcation diagram of this family we find coher- 
ence in the phase portraits when "traveling" from 
one part to the other. 

Analogously to the previous class, we shall foli- 
ate the 3-dimensional bifurcation diagram in R 3 by 
planes h = ho, ho constant. We shall give pictures 
of the resulting bifurcation diagram on these planar 
sections on an affine chart on R 2 . In order to de- 
tect the key values for this foliation, we must find 
the values of parameters where the surfaces have 
singularities and/or intersect to each other. We re- 
call that we will be only interested in non-negative 
values of h to construct the bifurcation diagram. 

The following set of eleven results study the 
singularities of each surface and the simultaneous 
intersection points of the bifurcation surfaces, or 
the points or curves where two bifurcation surfaces 
are tangent. 

We shall use the same set of colors for the bi- 
furcation surfaces as in the previous case. 

Lemma 7.2. Concerning the singularities of the 
surfaces, it follows that: 

(i) (Si) has a straight line of singularities on 
(0, Z, 1/2), I G R; 

(ii) (S2) has a straight line of singularities on 
(0,0,m), m G R; 

(Hi) (S3) and (S5) have no singularities beyond 
their multiplicity as double planes; 

(iv) (Sq) has a straight line of singularities on 
(0, Z, 1/2), / G R. 

Proof. The gradient vector of surface (Si) is given 
by VSi(h,l,m) = (8h(l + 3hl - 2m),8h 3 , -8h 2 ), 
and solving the equation VSi(h,l,m) = (0,0,0) 
we get that this surface has a straight line of sin- 
gularities on (0, Z, 1/2), / G R, proving (i). The 
surface (62) trivially has a straight line of singu- 
larities on (0,0, m), m G R, which proves (ii). As 
the surfaces (53) and (£5) are two double planes, 
they have no singularities; so (Hi) is proved. To 
prove (iv), we see that surface (Sq) is the product 
of a plane and a quartic, each one not having singu- 
larities itself. However, their intersection produces 
a straight line of singularities for this surface on 
(0, Z, 1/2), / G R. ■ 
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Lemma 7.3. If h = 0, the surfaces (Si) and (S2) 
coincide. For all li / 0, they intersect along the 
straight line (h, 0,1/2). 

Proof. By restricting both equations of surfaces 
(1S1) and (1S2) to h = they become both null, coin- 
ciding on the plane h = 0. For all h 7^ 0, by solving 
simultaneously both equations of these surfaces, we 
obtain the straight line (h, 0, 1/2). ■ 

Lemma 7.4. The surfaces (Si) and (S3) coincide 
on the plane h = and have no intersection for all 
h^0. 

Proof. By restricting both equations of surfaces 
(<Si) and (53) to h = they become both null, co- 
inciding on the plane h = 0, and it is easy to verify 
that both equations have no simultaneous solutions 
for all h 0. ■ 

Lemma 7.5. If h = 0, the surfaces (Si) and (S§) 
intersect along the straight line (0, Z, 1/2), I G R. 
For all h 7^ 0, they intersect along the straight line 
(h, 0,1/2). 

Proof. By restricting both equations of surfaces 
(Si) and (Sq) to h = and solving the restricted 
equations simultaneously we obtain the straight line 
(0, 1, 1/2), / G R, and by solving the system formed 
by the equations of these surfaces for all h ^ 0, 
we find the straight line (h, 0,1/2). Moreover, we 
see that both straight lines intersect at the point 
(0,0,1/2). ■ 

Lemma 7.6. If h = 0, the surfaces (Si) and (Sq) 
coincide and, if h^ 0, they have a 2-order contact 
along the straight line (h,0, 1/2). 

Proof. By applying the same technique as in 
Lemma 17.31 we have the result. To prove the 2- 
order contact of the surfaces along the straight line 
(h, 0,1/2), h ^ 0, we note that (Si) n (<S 6 ) = 
{(0,1, m), l,m G R} U {(h, 0,1/2), lei}. Then, 
the gradient vector of both surfaces along the 
straight line (h, 0,1/2) is a multiple of the vector 
(0, h, —1), /i / 0. If we check the next derivative, 
they do not coincide, which proves the 2-order con- 
tact. ■ 

Lemma 7.7. The surfaces (S2) and (S3) coincide 



on the plane h = and have no intersection for all 
h^0. 

Proof. Analogous to Lemma 17.41 ■ 

Lemma 7.8. The surfaces (S2) and (S§) inter- 
sect along the straight lines (0,1,1/2), I G R and 
(h, 0,1/2), h^0. 

Proof. Analogous to Lemma 17.51 ■ 

Lemma 7.9. The surfaces (S2) and (Sq) coincide 
on the plane h = and intersect along the straight 
line (h, 0,1/2), h^0. 

Proof. By restricting both equations of surfaces 
(52 ) and (Sq) to h = they become both null, co- 
inciding on the plane h = 0, and by solving the sys- 
tem formed by the equations of these surfaces for 
all h 7^ 0, we find the straight line (h, 0, 1/2). ■ 

Lemma 7.10. The surfaces (S3) and (S§) inter- 
sect along the straight line (0,1, 1/2), / G BL 

Proof. By solving the system formed by the equa- 
tions of these surfaces we find the straight line 
(0,1, 1/2), I e R. ■ 

Lemma 7.11. The surfaces (S3) and (Sq) coincide 
on the plane h = and have no intersection for all 
h^0. 

Proof. Analogous to Lemma 17.41 ■ 

Lemma 7.12. The surfaces (S§) and (Sq) inter- 
sect along the straight lines (0,1,1/2), I E R, and 
(h, 0,1/2), h G R, and the curve (-2/1,1,1/2), 
I ^ 0. 

Proof. By solving the system formed by the equa- 
tions of these surfaces we find the straight lines 
(0, Z, 1/2), I G R, and (h, 0,1/2), h G R, and the 
curve (-2/1,1, 1/2), 0. U 

Now, we shall study the bifurcation diagram 
having as reference the values of h where significant 
phenomena occur in the behavior of the bifurcation 
surfaces. As there is not any other critical value of 
h, except h = 0, this is the only value where the be- 
havior of the bifurcation surfaces changes critically. 
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Fig. 11. Slice of parameter space for © when h = 1. 



Fig. 10. Slice of parameter space for (J5} when h = 0. 

Recalling we are considering only non-negative val- 
ues of h, we shall choose a positive value to be a 
generic case. 

We take then the values: 



The value ho correspond to an explicit value of 
h for which there is a bifurcation in the behavior of 
the systems on the slices. The value h\ is just an 
intermediate point we call by a generic value of h 
(see Figs. [TUl and fTTj) . We recall that the bifurca- 
tion T = (two other finite points collide), plotted 
usually in green, does not appear here because it 
corresponds completely to the plane h = 0. 

As in the previous family, all the bifurcation 
surfaces intersect on h = 0. In fact, the equations 
of surfaces (Si), (S2), (S3) and (Sq ) are identically 
zero when restricted to the plane h = 0, and the 
equation of (1S5) is the straight line — 1 + 2m = 0, 
for all h > and m.Iel. 

Here we also give topologically equiv- 



alent figures to the exact drawings of 
the bifurcation curves. The reader may 
find the exact pictures in the web page 
http:/ /mat .uab.es /~artes /articles / qvfsn2SN02 / 
qvfsn2SN02.html. 

If we consider the value h = 1, other changes 
in the bifurcation diagram happen. On this plane, 
surface (Si) is the straight line 1 + 21 — 2m = 0, 
which intersects surface (1S5) at the point (1, 0, 1/2); 
surface (Sq) is the parabola 1 + 21 + I 2 — 2m = 
passing through the point (1,0, 1/2) with a 2-order 
contact with surface (<Si), see Lemma [Tl)l moreover, 
Lemma 17.121 assures that surface (Sq) has another 
intersection point with surface (£5) at (1,-2,1/2); 
surface (52) is the straight line / = 0, which in- 
tersects surfaces (<Si), (S$) and (Sq) at the point 
(1,0, 1/2), see Lemmas 17.81 and I7T91 and finally, sur- 
face (S3) is a negative constant. 

We recall that the black surface (56 ) (or W3) 
means the turning of a finite antisaddle from a node 
to a focus. Then, according to the general results 
about quadratic systems, we could have limit cycles 
around such focus for any set of parameters having 
W 3 < 0. 

In Figs. [T2] and [T3] we show the complete bi- 
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Fig. 12. Complete bifurcation diagram of QsnSN(B) 
for slice h = 0. 

furcation diagrams. In Sec. [9] the reader can look 
for the topological equivalences among the phase 
portraits appearing in the various parts and the se- 
lected notation for their representatives in Fig. [3l 

8. Other relevant facts about the bifurca- 
tion diagrams 

The bifurcation diagrams we have obtained for the 
families Q~snSN(A) and QsnSN(B) are completely 
coherent, i.e., in each family, by taking any two 
points in the parameter space and joining them by 
a continuous curve, along this curve the changes in 
phase portraits that occur when crossing the dif- 
ferent bifurcation surfaces we mention can be com- 
pletely explained. 

Nevertheless, we cannot be sure that these bi- 
furcation diagrams are the complete bifurcation di- 
agrams for Q~snSN(A) and Q~snSN(B) due to the 
possibility of "islands" inside the parts bordered by 
unmentioned bifurcation surfaces. In case they ex- 
ist, these "islands" would not mean any modifica- 
tion of the nature of the singular points. So, on the 
border of these "islands" we could only have bifur- 
cations due to saddle connections or multiple limit 
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Fig. 13. Complete bifurcation diagram of QsnSN(B) 
for slice h = 1. 



cycles. 

In case there were more bifurcation surfaces, 
we should still be able to join two representatives 
of any two parts of the 66 parts of Q~snSN(A) or 
the 30 parts of Q~snSN(B) found until now with 
a continuous curve either without crossing such bi- 
furcation surface or, in case the curve crosses it, it 
must do it an even number of times without tan- 
gencies, otherwise one must take into account the 
multiplicity of the tangency, so the total number 
must be even. This is why we call these potential 
bifurcation surfaces " islands" . 

However, in none of the two families we have 
found a different phase portrait which could fit 
in such an island. The existence of the invariant 
straight line avoids the existence of a double limit 
cycle which is the natural candidate for an island 
(recall the item (jivj) of Sec. H|), and also the limited 
number of separatrices (compared to a generic case) 
limits greatly the possibilities for phase portraits. 
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9. Completion of the proofs of Theorems 
[TTQandCLl 

In the bifurcation diagram we may have topolog- 
ically equivalent phase portraits belonging to dis- 
tinct parts of the parameter space. As here we 
have finitely many distinct parts of the parameter 
space, to help us identify or to distinguish phase 
portraits, we need to introduce some invariants and 
we actually choose integer-valued invariants. All 
of them were already used in |Llibre et al., 2004 
Artes et al., 2006] . These integer-valued invariants 
yield a classification which is easier to grasp. 

Definition 9.1. We denote by I\(S) the number 
of the real finite singular points. This invariant is 
also denoted by f% j(S) |Artes et al., 2006| . 

Definition 9.2. We denote by h(S) the sum 
of the indices of the real finite singular points. 
This invariant is also denoted by deg(DIf(S)) 
| Artes et al, 2006| . 

Definition 9.3. We denote by ^(S 1 ) the number of 
the real infinite singular points. This number can 
be oo in some cases. This invariant is also denoted 
byN Mi00 (S') | Artes et a/.,"2006l . 

Definition 9.4. We denote by I4(S) the sequence 
of digits (each one ranging from to 4) such that 
each digit describes the total number of global or 
local separatrices (different from the line of infinity) 
ending (or starting) at an infinite singular point. 
The number of digits in the sequences is 2 or 4 
according to the number of infinite singular points. 
We can start the sequence at anyone of the infinite 
singular points but all sequences must be listed in 
a same specific order either clockwise or counter- 
clockwise along the line of infinity. 

In our case we have used the clockwise sense 
beginning from the saddle-node at the origin of the 
local chart U± in the pictures shown in Figs. [I] and 
[21 and the origin of the local chart 17% in the pictures 
shown in Fig. [3j 

Definition 9.5. We denote by -^(S 1 ) a digit which 
gives the number of limit cycles. 

As we have noted previously in Remark 16.181 



we do not distinguish between phase portraits 
whose only difference is that in one we have a fi- 
nite node and in the other a focus. Both phase 
portraits are topologically equivalent and they can 
only be distinguished within the C 1 class. In case 
we may want to distinguish between them, a new 
invariant may easily be introduced. 

Definition 9.6. We denote by Iq(S) the digit 
or 1 to distinguish the phase portrait which has 
connection of separatrices outside the straight line 
{y = 0}; we use the digit for not having it and 1 
for having it. 

Definition 9.7. We denote by It(S) the sequence 
of digits (each one ranging from to 3) such that 
each digit describes the total number of global or 
local separatrices ending (or starting) at a finite 
antisaddle. 

Theorem 9.8. Consider the subfamily QsnSN(A) 
of all quadratic systems with a finite saddle-node 
sn(2) an d an infinite saddle-node of type (°) SN lo- 
cated in the direction defined by the eigenvector with 
null eigenvalue. Consider now all the phase por- 
traits that we have obtained for this family. The val- 
ues of the affine invariant I = (I±, I2, 13, 14, 15, Iq) 
given in the following diagram yield a partition of 
these phase portraits of the family QsnSN(A). 

Furthermore, for each value of I in this dia- 
gram there corresponds a single phase portrait; i.e. 
S and S' are such that I(S) = I(S'), if and only if 
S and S' are topologically equivalent. 

Theorem 9.9. Consider the subfamily 
QsnSN(B) of all quadratic systems with a fi- 
nite saddle-node sum) and an infinite saddle-node 

of type (2) SN located in the direction defined by 
the eigenvector with non-null eigenvalue. Consider 
now all the phase portraits that we have obtained 
for this family. The values of the affine invariant 
1 = I2, 13, h, It) given in the following diagram 
yield a partition of these phase portraits of the 
family QsnSN(B). 

Furthermore, for each value of I in this dia- 
gram there corresponds a single phase portrait; i.e. 
S and S' are such that I(S) = I(S'), if and only if 
S and S' are topologically equivalent. 
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The bifurcation diagram for Q~snSN(A) has 66 
parts which produce 29 topologically different phase 
portraits as described in Table 1. The remaining 
37 parts do not produce any new phase portrait 
which was not included in the 29 previous ones. 
The difference is basically the presence of a strong 
focus instead of a node and vice versa. 

Similarly, the bifurcation diagram for 
QsnSN(B) has 30 parts which produce 16 
topologically different phase portraits as described 
in Table 3. The remaining 14 parts do not produce 
any new phase portrait which was not included 
in the 16 previous ones. The phase portraits 
are basically different to each other under some 
algebro-geometric features related to the position 
of the orbits. 

The phase portraits having neither limit cy- 
cle nor graphic have been denoted surrounded by 
parenthesis, for example (Vi) (in Tables 1 and 3); 
the phase portraits having one limit cycle have been 
denoted surrounded by brackets, for example [Vn] 
(in Table 1); the phase portraits having one graphic 
have been denoted surrounded by {}, for example 
{751} (in Table 1). 

Proof. The above two results follow from the results 
in the previous sections and a careful analysis of the 
bifurcation diagrams given in Sees. 0and[71 in Figs. 
El O [12] and [131 the definition of the invariants Ij 
and their explicit values for the corresponding phase 
portraits. ■ 

We first make some observations regarding the 
equivalence relations used in this paper: the affine 
and time rescaling, C 1 and topological equivalences. 

The coarsest one among these three is the topo- 
logical equivalence and the finest is the affine equiv- 
alence. We can have two systems which are topo- 
logically equivalent but not C 1 -equivalent. For ex- 
ample, we could have a system with a finite anti- 
saddle which is a structurally stable node and in 
another system with a focus, the two systems being 
topologically equivalent but belonging to distinct 
C 1 -equivalence classes, separated by the surface Sq 
on which the node turns into a focus. 

In Table 2 (Table 3, respectively) we listed in 
the first column 29 parts (16 parts, respectively) 
with all the distinct phase portraits of Figs. [Hand 
M (Fig. respectively) . Corresponding to each part 
listed in column 1 we have in its horizontal block, all 



parts whose phase portraits are topologically equiv- 
alent to the phase portrait appearing in column 1 
of the same horizontal block. 

In the second column we have put all the parts 
whose systems yield topologically equivalent phase 
portraits to those in the first column, but which 
may have some algebro-geometric features related 
to the position of the orbits. 

In the third (respectively, fourth, and fifth) col- 
umn we list all parts whose phase portraits have 
another antisaddle which is a focus (respectively, a 
node which is at a bifurcation point producing foci 
close to the node in perturbations, a node-focus to 
shorten, and a finite weak singular point). In the 
sixth column of Table 1 we list all phase portraits 
which have a triple infinite singularity. 

Whenever phase portraits appear on a horizon- 
tal block in a specific column, the listing is done 
according to the decreasing dimension of the parts 
where they appear, always placing the lower dimen- 
sions on lower lines. 

9. 1 . Proof of the main theorem 

The bifurcation diagram described in Sec. El plus 
Table 1 of the geometrical invariants distinguish- 
ing the 29 phase portraits, plus Table 2 giving 
the equivalences with the remaining phase portraits 
lead to the proof of the main statement of Theorem 
11.11 Analogously, we have the proof of Theorem 
11.21 but considering the description in Sec. [7] and 
Tables 3 and 4. 

To prove statements (c) and (d) of Theorem 
11.21 we recall the Main Theorem of |Vulpe, 2011| 
and verify that: 

(i) Any representative of 2S\ is such that h > 0, 
I = and m > 1/2. Then, we have: Ta = 0, 
T 3 = 8h 2 (-l + 2m) ^ 0, T Z T = -8/i 4 (-l + 
2m) 3 < 0, T\ = Ti = = 0, which imply 
that it has an integrable center c; 

(ii) Any representative of 2S2 is such that h > 0, 
I = and m < 1/2. Then, we have: T4 = 0, 
T 3 = 8/i 2 (-l + 2m) + 0, %T = -8/i 4 (-l + 
2m) 3 > 0, F\ = Ti = = 0, which imply 
that it has an integrable saddle $. 
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Table 2: Topological equivalences for the subfamily QsnSN(A). 
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Table 3: Geometric classification for the subfamily QsnSN(B). 
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